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Transformations  of  Gaussian  Random  Fields  and  a 
Test  for  Independence  of  a  Survival  Time 
from  a  Covariate 

Ian  W.  McKeague  A.  M.  Nikabadze  Yanqiiig  Sun 


Abstract 

It  has  been  almost  sixty  years  since  Kolmogorov  introduced  a  distribution-free  om¬ 
nibus  test  for  the  simple  null  hypothesis  that  a  distribution  function  coincides  with  a 
given  distribution  function.  Doob  subsequently  observed  that  Kolmogorov’s  approach 
could  be  simplified  by  transforming  the  empirical  process  to  an  empirical  process  based 
on  uniform  random  variables.  Recent  use  of  more  sopliisticated  transformations  has  led 
to  the  construction  of  asymptotically  distribution-free  omnibus  tests  when  unknown 
parameters  are  present.  The  purpose  of  the  present  paper  is  to  use  the  transformation 
approach  to  construct  an  asymptotically  distribution-free  omnibus  test  for  indepen¬ 
dence  of  a  survival  time  from  a  covariate.  The  test  statistic  is  obtained  from  a  certain 
test  stJ.tistic  process  (indexed  by  time  and  covariate),  which  is  shown  to  converge  in 
distribution  to  a  Brownian  sheet.  A  simulation  study  is  carried  out  to  investigate  the 
finite  sample  properties  of  the  proposed  test  and  an  appUcation  to  data  from  the  British 
Medical  Research  Council’s  4th  myelomatosis  trial  is  given. 


1  Introduction 

A  standard  way  of  testing  for  independence  of  a  survival  time  from  a  covariate  z  is  to  fit 
Cox’s  (1972)  model  for  the  conditional  hazard  function,  A(flz)  =  Ao(<)e^^,  and  test  whether 
the  regression  parameter  ^  is  zero.  However,  this  test  has  limited  power  because  of  the 
restrictive  (viz  parametric  and  multiplicative)  modeling  of  the  covariate  effect. 

In  this  paper  we  develop  an  omnibus  test  that  can  detect  arbitrary  forms  of  dependence 
of  a  (possibly  censored)  survival  time  on  a  one-dimensional  covariate,  and  which  is  asymptot¬ 
ically  distribution-free.  The  latter  property  wdll  be  achieved  via  the  transformation  method 
of  Doob  (1949)  and  Khmaladze  (1981). 

We  begin  by  giving  some  background  to  the  general  problem  of  constructing  omnibus  tests 
(i.e.  tests  consistent  against  all  alternatives)  which  have  the  distribution-free  property.  First 
consider  the  simple  hypothesis  F  =  Fq,  where  Fo  is  specified  and  the  life  times  Tj,  •  •  • ,  Tn 
are  completely  observed  iid  random  variables  having  distribution  function  F.  Let  F{i)  = 

^AMS  1991  subject  dassi/iCHtlons.  Primary:  G2G 10;  secondary:  C2G20 

^Key  words  and  Phrases.  Distribution-free  omnibus  test,  innovation  Brownian  sheet,  counting  proccssc.s, 
martingale  methods. 
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fiT,  <  t)  be  tlie  empirical  distribution  function  of  tlie  7','s  and  e„(0  “  v'^(^'(0  '■ 
Fo(<))  tlie  eni])irical  process.  Assume  that  /'o  is  contiimoii.-i.  Doob  (1919)  liaiisformed  i’n(l) 
to  the  uniform  empirical  i)rtKes8  u,,(r)  —  'vlucb  is  an  empirical  prt>coss  based 

on  the  iid  uniform  random  variables  /‘o(V’i),  i  =  l.-'-.n.  The  distribution  of  ?/„  dot's  not 
depend  on  Fq  (and  it  converges  weakly  to  a  Brownian  bridge),  so  the  di.V.iiimtion  tif  any  test 
statistic  that  is  a  functional  of  is  free  from  /'o.  In  particular,  tlie  Kolmogorov  Smirnov 
statistic  supj.  |«7i(.r)|  and  the  Cramer  von  Mises  statistic  j  iil(x)(Lr  are  distribution-free. 

Next  consider  the  composite  null  hypothesis  F  —  I'o(-,0).  where  0  is  an  unknown  pa¬ 
rameter.  The  natural  extension  of  the  above  transformation.  »„(.r)  =  i\,(l\'^\.r.O)),  where 
//„(/)  =  s/ri[F{l)  —  Fo{tJ}))  is  the  paranu'tric  «'mi)irical  procc'ss  and  0  is  an  estimator  of  0,  is 
unforturnately  no  longer  di.stribution-free  or  even  asymptotically  distribution-fri'c  (Durbin. 
1973).  -As  a.  consecpieuce,  classical  .statistics  such  as  su|)j.  |(7„(.r)|  or  /  u^j{.v)  dr  have'  limit 
distributions  which  depend  on  Fa-  'riius  it  is  necessary  to  construct  a  more  sophisticated 
transformation  of  On  that  can  provide  the  basis  for  gooduess-of-fit  tests,  geiu-ralizing  what 
the  uniform  empirical  proce.ss  does  in  the  case  of  simple  hypothe.ses.  Khmaladze  (1981) 
introduced  martingale  methods  to  address  this  problem;  sec'  also  Nikabadze  (1987).  riic' 
parametric  empirical  process  O,,  converges  weakly  to  some  zero-mean  (laussiau  ])rocess  /z 
(Durbin,  1973),  so  Khtnaladw  first  transformerl  the  proee'j^  </  to  an  innovation  rnarlingalr, 
which  is  a  Gaussian  process  with  iudcpeiident  iiu  rements  and  covariance  function  Ai,0) 
and  which  preserves  the  information  in  /z.  Then  he  transformed  the  innovation  martingale 
to  a  standard  Brownian  motion  in.  Applying  the  transformation  u  »-♦  w  to  j7,i,  lesult.-'  iii 
a  test  process  that  converges  weakly  to  Brownian  motion.  This  leads  to  an  asymptotically 
distribution-free  omnibus  test. 

Chi-squared  tests  are  widely  used  for  goodness-of-fit  testing  and  for  testing  independence 
of  two  variables  in  a  contingency  table  analysis.  Thc'v  were  first  introduced  by  Pearson  (1900) 
for  .simple  hypolhesi*'  F  —  /'«.  The  rhi-scpiarc'd  .statistic  is  Ibniusi  l>y  dividing  pari  of  the* 
real  line  into  cells  and  comparing  the  observed  and  <'xi)ected  frequency  in  each  cell.  Fisher 
(1922,  1924)  extended  this  statistic  to  handle  the  pre.^ence  of  an  unknown  parameter  0  in 
Fo-  Chi-squared  Utsls  depend  on  an  arbitrary  choice  of  intervals  and  tin  y  only  use  grouped 
data.  Although  chi-squared  tests  are  easy  to  ])erform,  they  are  not  omnibus  (unless  the 
variables  are  discrete)  and  are  tyjrically  less  powerful  than  tc'sts  of  Kolmogorov  Smirnov  or 
Cramer-von  Mises  type,  which  use  all  the  information  in  the  data. 

In  survival  analysis,  one  is  rarely  able  to  observe  complete  life  histories.  Important 
examples  occur  with  right  ceim^rihg  and  left  '•.•uucatU>ii  (Kcidiiig  and  GiH.  PHit)).  These 
examples  fit  into  the  general  setting  of  Aalen's  (1978)  rmiltiplicativc  intensity  model  for 
counting  processes.  In  that  setting  it  is  natural  to  formulate  hypotheses  iii  terms  of  the 
haZaid  fauctioli  A(f)  oi  tlie  cumulative  hazaid  fuiicliuu  A{1)  =  /q  A(s)  (is.  rather  than  the 
distribution  function  F.  .Andersen  ct  al.  (1982)  studied  tests  of  the  simple  hypothesis  A  =  Ao 
in  terms  of  fuuclioi.als  of  ■  -^u);  whore  A  is  the  Nelson -.A aUna  estimator,  lljort  (199J) 
considered  the  composite  hypothesis  A  =  Aof-.fD-  with  statistics  leased  on  functionals  of  the 
process  y/n{  A{f)  —  .'\o(i,^)),  where  0  is  the  maximum  likelihood  estimator  of  0.  'I'liis  process 
coioeiges  iseakis  tu  a  aciu-iUeali  Gaussian  proccs?.  uiidei  the  null  Inpolhesis,  and  can  be 
used  to  construct  chi-squared  tests.  Alternatively,  an  innovation  martingah'  can  be  found 
for  tlie  limit  imxh-oss  .™nd  used  consiruel  an  :»iyn!ptot.M-(il!y  distrilmthni  fro'  omnibus 

In  man\'  applications  of  survival  analysis  it  is  important  to  consider  wlicthei-  a  covariate 


has  some  effect  upon  survival,  say  through  the  conclitioual  hazard  luiictioii  A(/|;)  =  A(/.c). 
That  is.  one  would  like  to  test  the  mill  hypothesis 

//o  :  A(/,r)  does  not  depend  oil  the  (ovariate 

against  the  general  alternative  that  A(/,c)  deiieiids  on  c.  I'or  simplicity.  w('  shall  restrict 
the  domain  of  (/,;)  to  be  the  unit  srpiare.  An  oiimibiis  test  of  //(,  is  leasihh'  wlu'ii  the 
covariate  is  one-dimeiisional,  such  as  ago  at  diagnosis,  disease*  diiratioii,  etc.  Indeed,  .McK- 
eague  and  Utikal  (1990,  subsequently  MU)  jn-oposed  such  a  te'st  liast'el  on  the  process 
A'(/,-)  =  —  -4),  where  A  is  an  estimate  of  the  doubly  cuiiiiilat ive  hazai'l  function 

A{(,:)  =  JiJq  \(s.  x)  (I.T  (l.'i,  and  A{t.z)  =  -A(/)  is  the  natural  ('st  iiiiate  of  A  under  //d.  They 
showed  that  A’  converges  weakly  under  //(»  to  a  (iaiissiaii  random  !i<*l(l  of  the*  loriii 

m(/,c)=  /  /  \/Ti(l\\'  —  h{z)  f  f  <j  <I\V.  (1-1) 

JoJo  Jo.lo 

where  lU  is  a  nrowiiian  sheet  and  hjKfJ  are  certain  nonrandom  fuuctions.  The  aliovi*  stochas¬ 
tic  integrals  are  defined  in  the  L'^-sense,  see  Wong  and  Zakai  (1!)71).  Ml''s  test  was  leased  on 
the  Kolmogorov  Sirvinu-v  statistic  d  direct.ly  Ctom  A^  Honw^vt.  while  awmptotirally 

omnibus,  such  a  test  is  not  asymptotically  distribution-free  and  would  r<'(nure  simulation  of 
the  process  m  to  find  critical  values. 

We  siiall  construct  a  transfonnatioiv  J  that  iuai>s  m  t,<>  its  innovation  Urosvnian  shwt.  An 
estimated  version  J  of  J  will  be  obtained  by  plugging  an  estimate  of  li  int  o  .7  (it  turns  out  that 
J  does  not  involve  g  and  b  is  known).  We  then  show  that  J{X)  converges  weakly  to  Urownian 
sheet.  In  this  way  we  obtain  an  asymptotically  distribution-free  omnibus  test  for  JIq,  with 
the  Kolmogoiov-Smirnov  statistic  computed  from  J(X).  No  .sinnilntioii  techiiiciue  is  needed 
to  find  critical  values.  The  test  statistic  coiivru’ges  weakly  to  su])  \\V(I.  ;)|.  Alt  hough  an  exact 
formula  for  the  distribution  furicticii  of  sup  jlTf/.  .r)|  is  not  known  (only  apjrroximatious  are 
available,  see  Adler  (1991)),  it  is  straightforward  to  carry  out  a  single  Monte  Carlo  exireriment 
to  evaluate  it  quite  accurately.  Thus,  our  test  avoids  difficulties  arising  from  simulating  the 
null  distribution  for  each  particular  problem. 

The  paper  is  organized  as  follows.  In  Section  2.  we  construct  the  tiausformation  J.  In 
Section  3.  we  introduce  the  estimate  J  and  define  the  test  statistic.  Results  of  a  simulation 
study  are  reported  in  Section  -1.  In  Section  o,  the  test  is  ap])lied  to  a  set  of  data  from  the 
Hritish  Medical  Research  Council’s  (1981)  -Uli  myelomatosis  t  rial,  l^roirerties  of  the  test  are 
proved  in  Section  G.  \  arious  lemmas  needed  through  the  ])a])('r  are  collected  in  an  appendix. 


2  Transformation  of  m  to  Brownian  Sheet 

In  this  section  we  construct  our  transformation  ./  of  the  Gaussian  random  field  ui  in  (1.1)  to 
Hrownian  sheet.  Such  a  <  raiisformatioii  is  likely  to  have  further  applications  in  nonparametric 
statistics  beyond  uur  test  for  invlepeudeiur  •  in  .any  setting  where  a  test  process  vv.'nNa.'rgr.5 
weakly  to  a  process  of  the  form  (1.1):  e.g.  in  testing  whetlier  \{1.z)  is  independent  of  /  (i.e. 
the  roles  of  i  and  r  are  reversed),  or  testing  whether  a  jiure  jumi)  ])rocess  on  a  finite  state* 
Sjjtui  IS  rt  si-liii-Mdi  ken  ))Iocess,  sea'  Mb  ^SeC^lOll  1.2).  Df  eoiUse.  it  is  usii<ill\  iieie’sscU}  to 
estimate  ./  and  how  that  is  done  will  deijcnd  on  the  particular  application. 


We  begin  with  a  key  proposition  showing  that  the  law  of  a  Brownian  slu'('t  B'  is  prc'served 
under  a  shift  of  B’  by  a  eert.ain  fnnetional  of  B’. 

Proposition  2.1  Let  Ic  <=  //*’([(),  Ij-)  satisfy  k-(s.v]elr  >  t)  a.t  .  [c/s]  for  n  <  1 .  and  let  B' 
be  a  lirenvtrian  sheet.  Tliiu 


n{t,=)  =  \vm=)- 


'  r'  k{s.  .r)k{s.  ii) 


./()  l./o./u  f]  k'"{s.v)  elv 
is  a  Brownian  slixet  on  [0,1  ]‘.  The  relation  (2.2)  is  inn  rtilile : 

f'i‘‘  k{s..r)k(s.u) 


<l\V{s..r) 


(hi 


{•2.2) 


Proof  Let, 


Then 


ir(/.~')  =  /^(/.-')  + 


(l{t,V\S,.T)  = 


B{t.,z)  =  WiTz)- 


./o  L./o./„  //F(s.r)de 

k{.'^.x)k{s.n)I(.r  >  ii)I[s  <  t) 
/jF(.s.r)Jr 

1  /•! 


(In. 


a{t.  ir.s,  x)  dB’(.s.  .r) 


(In. 


Jo  L./0./0 

Notice  that  B  is  a  Gaussian  random  field,  so  we  only  need  to  inspect  its  covariance  funct  ion. 

For(^^^')€[0,l]^ 


cov(B(t,z}.B(f'.z'})=(tAf')(zAz') 

~  /  Iff  (i(t\ii':s,.r)(ls(l.r  elu'  —  f  \  f  f  (i{t  ,ir.  s..r)  ds  dr 
Jo  L  JqJo  J  ./()  L  Jo  Jo 


da 


-'/■j'  r  /■!  /■! 


+ 

JoJo  L./0./0 
=  (/A/')(-Ac') 


a{t.  a:  s.  x)a(t'.,  u'\  s.  .r)  (/■■.  d.r 


(III  dll' 


1  /-I 


+ 


.10.(0  LJo.'o 


a{l.u;  s.x)a{t\  ii\  s,  ,f )  ds  dx 


—  f  a(t.ir,s.n')l(.s  <  l')d.s  —  f  a(t\n'\s,ti)l{s  <  t)ds 
Jo  Jo 


dll'  (In. 


Since 


/J  «(/,  ir.  s.  x)a{l'.  ii'\  s,x)  dx 

k{s,ii)l{s  <  l)k(s,v')l{s  <  t')  f^  k‘(.s.x)/(.r  >  ii  V  ii')d:r 
!^,kHs.v)dvCkH.s.v)dv 
k{s,  u')k{s.ii)I{n'  >  ii)I{s  <  /)/(.s  <  /') 


+ 


k%‘(,  v)dv 

/:(.s,  u)/(ii  >  n')/(.s  <  t)l{s  <  I') 


Ck%^.v)dr 
=  a(l.  ir.s.  ii')l{s  <  I')  +  n{l'.  ti'.  s.ti)I{s  <  I). 


for  almost  all  (n.i/'.  >.)  ^  [0,1]’’,  we  have  that  B  is  a  Brownian  sheet.  It  can  be  verilied 
immediately  that  (2.3)  is  the  inverse  of  (2.2).  □ 


A  Brownian  motion  »>(/)  is  calK'd  an  iininvation  process  of  a  process  ^(/)  if  ir  carries 
the  same  “information''  as  the  ])rocess  i.e.  the  ^r-fields  JT,"'  and  generated  l)y  (<•  and  ^ 
up  to  each  time  /  coincide,  see  Liptscr  and  Shiryayev  (1977,  ]).  200).  For  onr  ]nn  poses.  the 
appropriate  extension  of  this  definition  to  a  two-i)arameter  process  and  a  Brownian 

sheet  is  made  by  requiring  where  =  (t{^(s,c);  c  £  [0.  l],.s  t  [h./]}  and 

fP  is  similarly  defined.  Note  that  rcpre.sents  the  information  about  at  all  values 

of  c  and  all  s  <i. 

We  now  give  the  main  result  of  this  section,  providing  an  innovation  Brownian  sheet  for 
the  ])rocess  in  (1.1). 


Theorem  2.1  Svpposc  ihai  h  :  [0.  1]‘  — >  Si  is  a  hovtulcd  posillrr  iii((istir(il)l(  fuiirlinii  which 
is  bounded  aicny  from  zero,  h  :  [0,  1]  — *  Si  is  dijfi  rcnliahlr  with  sijuarf  inlnirahh  drvivoiivr. 
fl{b'{x)pdx  >  b,  c  €  [0,  1 )  and  g  £  L‘([Q.  1]^).'  Then 


B{t,z)=  (  f  Ir^din—  [  f  f  Ir ^[s.v)Q{s,u.x)  dx 
JoJo  ./lUO  L./o 

is  an  intiovaiioii  Brownian  sheet  of  the  process  m,  where 


din{s.ri)  ('--‘1) 


Q{s,u,x)  = 


h  7  ( .s .  j/ )  //( H )  /i  -  ;  ( .s ,  :r )  /)' ( .j- ) 

Ir  h-^{s,v){b'{v)ydv 


Proof  Notice  that 

rl  ruAz 


Let 


f  [  Q{s,u,x)ir^{s.u)b'{n)dxdti  =  f  h~'^{s.x)h'(x)dx. 
Jo  Jo  Jo 

(/(/.)=  ffgdW. 

JoJo 


Substituting  m  into  (2.4)  we  get 


rt  rl  ruA: 

/(l./o  L  .'0 


+ 


Q{s,v,x)  dx 

ri'  \  r 

oJo  L  Jo 


f/li’(s.  (/) 

»)//((/)  dn  I  '{ils) 
d\\'{s,y} 


Qis.  11.  x)  dx 

=  ir(L.r)-  [[  g{s,y)  f  lr'^{s,x)l>{x)  dx 
JoJo  Ljo 

~  L  [lol 

+  11  [  /  Q(s.n.x)lr'^{s.v)b'(ii)dxdii  d\\'{s.y) 

JoJo  L  JoJo 


ir(/.c)- 


I  rl 


./()  L  i0.h- 


Q{s.  </,./•)  dll  ’(s.  //) 


dx. 


(•2.0) 


This  is  a  Brownian  sliec't  l)y  Lemma  2.1  with  h(s.x)  =  h  ^(.s.. /■)//(.;■).  The  innovation 
assertion  can  l)e  easily  obtained  by  (2. 4).  (2.5)  and  ap])Iication  of  Proposition  2.1.  □ 


We  shall  use  (ho  notation  ./  for  the  transformation  ^  i — ►  -/(O-  "'here  ^  is  a  random  lield 
and  ./(^)  is  defined  hy  the  right  side  of  (2.1)  with  to  re])lac('d  hy  I'lie  domain  of  ,/  is 
composed  of  random  fields  (  for  which  the  stochastic  integrals  in  ,/{(}  exist  in  the  A'-.sense. 
Theorem  2.1  shows  that  J{tn)  is  a  Hrowiiian  sheet. 


3  The  Test  Procedure 

In  this  section  we  first  desrrihe  the  roiinting  process  framework  h)r  onr  prolik'ni  and  fonnally 
define  A  and  A.  Then  we  show  that  the  transformation  ,/  given  al)ove  asymptotically 
transforms  A’  =  \fn{A  —  >1)  to  a  Mrownian  sheet.  'I'liis  is  doin'  \  ia  the  coiilinnons  mapping 
theorem,  h'inally,  we  construct  an  estimate./  ot  ./  and  show  that  J(\)  comi'iges  wi'akly  to 
a  Brownian  sheet.  This  will  complete  the  constrnclion  of  oni  test. 


3.1  The  Estimators  .4  and  A 

het.  N(/)  =  (A’l (/),•••,  jV,i(/)),  /  €  ^  multivariate  counting  procc’ss  with  res])ect  to 

a  right-continuous  filtration  i.c..  N  is  adapted  to  the  filtration  and  has  com])onents  N, 
which  arc  right-continuous  step  functions,  zero  at  timi'z<'ro,  with  jumps  of  size  -t-1  such  that 
no  two  components  jump  simultaneously.  .Assume  that  .W  has  intensity 


A.(/)  =  v;(/)A(/,z,(/)), 


where  is  a  predictable  {0,  l}-valued  process,  indicating  that  the  ftli  individual  is  at  risk 
when  )j(/.)  =  1,  and  Z,  is  a  predictable  [0,  l]-vahted  covariate  process.  The  function  A(/,.c) 
represents  the  failure  rate  for  an  individual  at  time  t  with  coviiriate  =  ~.  We  assume 
throughout  that  (A',-,  V’,.  Z,).  /'  =  1,  •  •  • . rt  are  iid  replicates  of  an  underlying  ' riph'  (.V.  ]  \  Z). 
Let  F{s^x)  =  P{Z3  <  .r,  Vj  =  1),  and  assume  that  for  each  s  g  [0,1].  is  absolutely 

continuous  on  [0,1]  with  suhdeiisity  The  functions  b.  /t,  (/  in  (1.1)  arc'  givc'ii  by 

b{z)  =  z,  li  =  A//  and  g  =  y/ A  ■  /.  The  transforniation  J  will  only  be  used  with  the.se  b  and 
h  from  now  on.  We  assume  that  /  and  A  are  Lipschitz.  of  houiidt'd  Viuiatioii,  and  bounded 
away  from  zero. 

Consider  t/„  equal  width  covariate  strata  —  [jr-i.-rv).  c  =  1. . .  .  wliere  .r,.  =  cu;,,. 
and  li'n  —  l/dn  is  the  stratum  width,  and  let  J;  =  for  r  G  I,  -  Tlic'u.  as  in  ML,  define 


A{Lz)  =  ff 


.V0‘>(c/s,.r 

V’0d(,s,  x) 


dx . 


(T6) 


where  N^"\t,z)  =  Jo  /(Z,(>>)  €  is  the  number  of  c-specilic  failures  olrscrved 

up  to  time  /,  and  Y'^’^\i.z)  =  l(Zi{t)  G  I:)Yi{i)  is  the  size  of  the  c-specifie  risk  .set  at 
time  t.  The  estimator  A  does  not  involve  stratificat  ion  of  the  covariate  and  can  bc'  obtained 
by  setting  ic„  =  1  in  A.  In  (3.G)  and  throughout  the  pajrer.  we  use  tin'  convc'iition  1/0  =  0. 


3.2  A  Continuous  Version  of  J 

\Vc  now  introduce  a  version  J  of ./  that  is  defined  on  a  snitalrly  large  function  sjrace  and  is 
continuous  on  a.  subspace  supporting  in.  so  the  coiitinuoiis  ma|)ping  theorem  is  applit  aide. 


(i 


Let  D'2  =  D2([0,  1]^)  l>o  llio  (’Xtonsioii  of  tlic  usual  Skoroliod  spaco  to  fuiuiious  oii 
[0,1]^.  see  Neuliaus  (1!)71).  I.e(  denote  tlie  siihspace  of  fuiiclions  ^  ^  l)>  fur  wli'xli 
have  hounded  variation,  and  let  (’2  <l<*note  lli<’  space  of  continuous  functions 
on  [0,1]^.  Kquip  1)2  and  ('2  with  the  uniform  norm. 

For  ^  €  C2  U  B\  2  i'nd  {t.~)  ^  [U.  I)  X  j0,/»).  with  ()</)<  1.  deliiu’ 


J(i){i.z)=  1 1  W^(.s..r )  -  1 1  (3.7) 

./0.'0  JoJo 

wdiere  the  integrals  are  coiisi<lered  to  he  w<"al<  integrals  (llildehraudt .  1!)()3).  ajul 


/2(.s.  «.r) 


h  Hs..r). 
/r'(.s.»)  / 


Jo 


fj  i')dv 


il.r . 


The  upper  hound  p  on  the  domain  of  r  is  used  to  ke(‘])  the  tlu'  <lenominator  in  f-}  hounded 
away  from  zero.  In  practice  p  would  h<*  taken  clos<'  to  1  (th<'  <'nd  of  the  range  <d  covariate 
values).  Note  that  J  is  a  w^ell-defined  map  from  C2U  Ii\  2  into  /.)2([0. 1]  x  [O./j])  since  Lemma 
1  shows  that  h  inherits  the  properties  of  /.A;  and  Lemma  2  ensures  the  e.\isteiice  of  the 
weak  integrals  when  ^  €  C'2.  We  have  include<l  IJ\2  in  th<*  domain  of  J  because  the  paths  of 
A'  belong  to  j3V^2i  hut  not  to  C2. 


Theorem  3.1  Suppose  (hot  w„  x  n~“.  where  f  <  o  <  1.  Then,  under  //q.  J{X)  com'crges 
weakly  to  a  Brownian  sheet  in  /52([0"  l]  x  [0, />)). 


Proof  Properties  of  /i,/2  obtained  via  Lemma  2  can  be  used  to  show  that  ./  is  continuous 
as  a  map  from  C’2  into  D2([0.  l]  x  [0,p]J.  In  particular,  we  use  the  pro])erty  that 
has  hounded  variation  uniformly  in  c,  0  <  ;  <  p.  MU  (Theorem  LI)  gives  that  A  converges 
weakly  in  D2  to  in.  where  in  is  defined  hy  (1.1)  with  b{z)  =  z.  Thus,  since  the  sample  paths  of 
m  belong  to  ('2  a.s..  the  continuous  mapping  theorem  (Billingsley.  lOfiSl)  gives  ./(A') — >.J{ni) 
in  Z)2([0. 1]  x  [O.p]).  The  processes  J{in)  and  J{in)  have  continuous  sample  paths  and.  by 
Lemma  3,  they  agree  a.s.  at  each  fi.xed  {t.z),  so  they  are  indistinguishable.  Theorem  2.1 
(with  b{z)  =  z)  implies  that  J{iu].  and  hence  J{tn).  is  a  Brownian  sh'*et.  □ 


3.3  Estimating  the  Transformation 

In  order  to  use  the  above  result  to  build  a  te.st  statistic,  we  need  to  estimate  the  unknowm 
function  in  J.  namely  li.  First  consider  the  kernel  estimator  h  suggested  by  .MU: 


h(i..z) 


t  -s 


z  —  x 


dn(s,x). 


w'here  is  .jandwidth  parameter.  A'  is  a  Lipschitz  nonnegative  ktMiu'l  function  with  com¬ 
pact  support  and  integral  1.  and 


/'/(/.  c  )  =  mi 


.V(")(J....r)  , 

- nr 


I 


is  an  estimator  of  //(/,  c)  =  f/JJ  /)(.s..r)f/sr/.r. 

\Vc  will  need  to  apply  methods  from  storliastic  ealciiliis  to  various  martingale  integrals 
involving  h,  which  is  possible  provided  that  is  an  -predict able  process  for  each  fixed 

c.  Since  />(•--)  is  continnoiis.  it  is  enough  that  it  be  adaptei'  to  the  filtration  jT,.  I'lnis.  we 
shall  use  a  kernel  function  I\  having  nonnegative  (as  well  as  compact )  support. 

The  estimated  transformation  ./  is  rlefined  by  inserting  a  truncated  version  h  of  h  in  place 
of  h  in  J,  where  h  is  given  by 

l,{L=)  =  h{l.=)l{c:'  <h(l.z)<c„). 

c„  >  0.  Note  that  ,7(A  )  is  well-defined  since  the  paths  o(  .\  belong  to  H\ 

3.4  The  Test  Statistic 

If  we  show  that  .7(A  )  converges  weakly  to  a  lirownian  shec't.  tlu-n  onr  test  for  //(,  can  be 
based  on  the  Kolniogorov-Smirnov  statistic 

.S’  =  snpo<,<,  M(\/»(i  -  >i))(f.  ;)1 

with  P-valnes  calculated  from  the  distribution  of  .snp(, ]ll'(/.c)].  For  that  purpose 
we  restrict  the  choice  of  u)„,6„,c„  as  follows: 

Condition  3.1 

V'n  X  71  where  -^  <  o  <  1. 
h„  X  7;“'’  where  0  <  J  <  |(1  -  o). 
c„  X  (log 77 1"*  where  *)  >  0. 

This  condition  is  satisfied,  for  example,  by  xi:„  x  h„  x  itK 

Theorem  3.2  Undrv  IIq,  J(X)  convcrgcf^  wrakly  lo  n  Hrmriiinii  slnct  in  /)j([().  1]  x  [O.p]), 

Onr  final  result  shows  that  the  lest  based  on  >’  is  oiiiiiibns.  consistent  agjiinst  any  depar¬ 
ture  from  the  null  hypothesis  Hq. 

Theorem  3.3  The  test  based  on  S  is  consistevt  against  the  general  alternatire  that  A(/.c) 
depends  on  z.  for  ij.z)  in  the  domain  [().  1]  x  [O.p]. 

4  A  Simulation  Study 

W'e  have  carried  out  a  limited  simulation  study  to  assess  the  performance  of  the  projroscd 
test.  We  considered  the  Kolmogorov-Sminiov  .statistic  .S'  with  the  supreiiumi  talo'ii  over 
[0,1]  X  [0.  .9].  i.e.  p  —  .9.  The  covariate  was  tak«'n  to  be  uniformly  distributed  ow-r  [0.1]. 
The  censoring  was  sinii)le  right,  censoring,  independent  of  the  failure'  time.  ;iiid  exironeiit ially 
distributed  with  parameter  adjusted  to  give  a  prescribed  irerceiitage  X/v  (low),  10%  (moder¬ 
ate),  and  75%  (heavy)  of  censored  observations  before  the  end  of  follow-nj).  Tin'  covariate 
strata  were  arranged  to  contain  etiual  numbers  of  observations.  For  sami)l('  size-  u  —  180, 


Tabic  1:  Observed  levels  and  powers  of  tlie  test  for  iiidc'jH'iideiiee  of  a  survival  time  from  a 
covariate  with  nominal  lev<>l  ’/a. 


A(/,r) 

samph' 

size 

censoring 

5% 

40'a 

75% 

ISO 

.038 

.021 

.0(10 

1 

300 

.012 

.039 

.003 

500 

.073 

.0.58 

.001 

ISO 

.008 

.025 

.007 

r‘ 

300 

.181 

.110 

.011 

500 

..525 

.38 

.009 

180 

.128 

.053 

.008 

300 

.118 

.279 

.032 

500 

1.872 

.708 

.182 

300  and  500,  the  mimber  of  strata  d„  was  taken  to  be  12.  15  and  20.  resulting  in  15.  20.  25 
covariate  values  in  each  stratum.  'I’lu*  corresponding  bandwidtlis  h„  W('n'  .20.  .20.  .22.  and 
the  kernel  function  A'  was  taken  to  be  the  indicator  of  [0. 1]. 

The  survival  times  were  generated  using  the  (ox  model  A(/,;)  =  for  Jo  =  0  (null 
hypothesis),  1  and  2  (alternative  hypotheses).  Table  4  gives  observed  levels  and  powers  of 
the  test  at  a  nominal  (a,sytnptotic)  level  5%,  with  each  entry  based  on  1000  samples.  In 
order  to  obtain  the  asymptotic  5%  critical  level  for  our  test  (i.e..  the  05th  p('rcentile  of 
®'^'Po<Ki.o<:<.9  1^^  (A  ')!)•  generated  10.000  replicates  of  the  Brownian  sheet  evaluated  on 
a  grid  defined  by  300  equally  spaced  points  on  each  axis,  idiots  of  the  density  and  disi  ribution 
function  of  the  supremum  of  the  absolute  value  of  the  Brownian  sheet  over  (0. 1]  x  [0.  .9]  are 
given  in  Figure  1.  The  5%  critical  level  was  found  to  be  2.28. 

The  simulation  results  show  that  the  observed  levels  arc  close  to  their  nominal  5%  values 
when  the  sample  size  is  at  least  ISO  and  censoring  is  light  or  moderate.  1  he  power  n'aches 
.52%  at  sample  size  500  and  low  censoring,  when  the  alternative  is  A(/.:)  =  (^  It  exceeds 
70%  at  sample  size  500  when  the  alternative  is  A(/.r)  =  and  cen.soriiig  is  inoch'rate. 

In  Figures  2-4  we  give  plots  of  the  observed  densities  (each  based  on  1000  sami)les)  of 
the  Kolomogorov-Sinirnov  statistic  S  under  the  null  and  alternative  hypotheses,  l  liey  are 
compared  with  the  density  of  snpy<.,<.m<;.<.  y  |ir(/.  c)j.  When  the  sample  size  is  at  least  300 
and  the  censoring  is  light  or  moderate,  the  observed  densities  agree  well  with  tlu'ir  theoretical 
limit  (see  Figure  2).  Under  the  alternatives  A(/,  r)  =  c'  and  c^',  wdien  the  sample  size  is  at 
least  500  and  the  censoring  is  light  or  moderate,  the  two  curves  are  quite  separate,  giving 
some  idea  of  the  power  of  the  test  (see  Figures  3  and  4).  In  I'igure  5.  we  give  perspective 
plots  of  a  realization  of  Brownian  sheet  and  a  realization  of  the  test  j)rocess  ./(.V)  (with 
A(Ac)  =  1.  sample  size  500  and  light  censoring).  As  ex])ected.  these  plots  are  qualitatively 
very  similar  to  one  another. 
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00 
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12  3  4 


(b) 


Figure  1:  (a)  The  clcii.sitv  and  (b)  the  distribution  function  of  sup  |lh’(f.c:)|,  {t,z)  E  [0, 1]  x 
[0,.9]. 
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Figure  2:  Observed  density  (dotted  line)  of  the  Kolmogorov-Smirnov  statistic  5  under  the 
null  hypothesis  X{t,:)  =  1  compared  witli  the  density  (solid  line)  of  sup  |lV'(<,c)l. 
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sample  size=500  sample  sjze=300  sample  size=1 80 


Figure  3:  Same  as  Figure  2,  but  for  the  alternative  hypothesis  A(t,c)  =  e". 
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sample  size=500  sample  size=300  sample  size=1 80 


1 3 


Figure  5:  A  single  rcjilization  of  Brownian  sheet  (first  row),  and  a  single  realization  of  J{X) 
(second  row)  from  simulated  data  {n  =  500,  light  censoring)  under  the  null  hypothesis 
A(t,2)  =  1.  Positive  parts  are  on  the  left,  negative  parts  on  the  right. 


\\ 


Figure  6;  The  lest  process  J(X)  lor  the  BMRC  data  (jrositive  part  on  the  lelt,  negative  part 
on  the  right). 
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5  Application  to  Myelomatosis  Data 

We  applied  our  test  to  a  set  of  data  from  the  British  Medical  Hesearcli  Couiicirs  (HMlKi) 
(1984)  4th  myelomatosis  trial.  The  data  set  contains  records  for  49.!)  j)atients,  iiicluiliug 
censoring  indicator,  scrum  microglobulin  (at  presentation)  and  survival  time  (in  days). 

Many  studies  (e.g.,  Cuzick  et  al.  (1985))  have  stiggested  that  serum  /?2  microglol)urm  has 
a  strong  effect  on  survival,  at  least  in  the  first  two  years  of  follow-up.  In  our  analysis  of  the 
data,  we  shall  ignore  all  covariates  except  for  serum  microglobulin,  which  is  taken  on  a 
login  scale  normalized  to  the  interval  (0,1).  The  end  of  follow-up  is  taken  to  he  2000  days, 
before  which  3%  of  the  observations  are  censored.  81  patients  were  still  at  risk  at  the  end 
of  follow-up.  The  survival  time  is  divided  by  2000  to  normalize  it  to  the  interval  [0, 1].  The 
covariate  interval  (0, 1]  is  divided  so  that  each  stratum  contains  20  covariate  values  except 
for  the  last  stratum.  We  used  p  =  .9,  as  in  the  simulation  study. 

We  have  plotted  the  test  process  J{X)  in  Figure  6.  The  magnitude  of  the  negative  part 
of  J{X)  shows  strong  departure  from  a  Brownian  sheet,  cf.  Figure  b.  'I'he  statistic  5‘  was 
found  to  be  5.335,  which  is  highly  significant.  Our  analysis  confirms  the  strong  influence  of 
a  patient’s  serum  ^2  microglobulin  on  survival. 


6  Proofs 


In  this  section  we  prove  Theorems  3.2  and  3.3.  We  begin  by  introducing  some  notation.  Let 
Mi  denote  the  .Frmartingale  A/,(<)  =  iV,(<)  -  Jq  A,(s)  ds,  and  set 


a(")(<,2)  =  x:a^.(^)€I.)i;(oa(<,^,(0)- 


<=i 


For  a  process  ^{t,  z),  set  ^^(0  =  ^r)  where  ir  =  ^Wn,  r  =  1 , . . . , 

We  shall  have  frequent  use  for  the  following  bounds  from  MU  (Lemma  1): 


sup  E 

5,x,n 


71W„ 


<  oc,  for  any  positive  integer  A’, 


sup /’(}'’*’‘^(s,  x)  =  0)  <  for  some  C  >  0. 


(6.S) 

(6.9) 


Proof  of  Theorem  3.2  By  Theorem  3.1,  it  is  sufiicienl  to  show  that  under  Ho, 

||(J- j)(A')||-t^0, 


where  II  •  II  is  the  supremum  norm  on  Z)2([0, 1]  x  [0,p]).  This  will  be  done  in  the  following 
two  steps: 


(G.IO) 

(G.ll) 


IG 


where  fi  =  /,  —  /,,  and  /■  is  obtained  by  inserting  h  in  place  of  h  in  /,.  i  =  1.2. 
Step  1  By  the  decomposition  of  X  given  in  MU  (proof  of  Theorem  -1.1 ), 

fj'hdx  = 

A(”)(5,x) 


(('*.12) 


— vn 


0  Jo 


+ /X  "0 

Ms,x)X(s)liY^''\s)  =  0)(ixds. 


where  M("l,  V’("l  are  defined  by  setting  =  [0, 1]  in  A/*"',  V*"',  resi)cctively.  VVe  denote  the 
four  terms  in  the  above  decomposition  by  /i,  /j,  /a,  /<,  respectively.  .Since  /\  is  contimious 
and  has  nonnegative  support,  we  have  that  ft(-,x),  and  therefore  /i(-,x),  is  J’rprcdictahle. 
Thus  the  stochastic  integrals  involved  in  Ii  and  are  square  integrable  martingales.  Now 
fill  is  bounded  by 

r*  rt  M(")(ds) 


where 


T){t)  =  sup 

><><<'» 


sup77(t)  +  v/nsup  I  /  /  /i  "~,V)7~f 


J  \  A/'"h 

£^(t,r)  and 


(6.13) 


is) 


Since  T;(t)  is  a  positive  submartingale,  Doob’s  inequality  gives  Esup,7;^(t)  <  4£’7/^(l).  Also, 
since  £^(l,r)  =  0,  and  E^(l,;)^(l, fc)  =  0  for  all  I  <  j  ^  k  <  dn,  we  can  apply  NlenchofF’s 
inequality  (see,  e.g.,  Shorack  and  VVellncr,  1986)  here  to  get 


<  0{\ogd„ftJ^  E[[Jj,(s.x)d:r) 

<  0(logd„)^^/'£[/  /i^(s,x)da:-^]ds 

JO  Ulr  Jr  (5)J 


ds 


<  0(logd„)V|  Jj 


=  0(logd,.)'[  /  /  E\frfdsdx\\ 

I  JO  JO 

The  second  term  in  (6.13)  is  hounded  by 

'■  ’'‘"/IMl  ^'-WTs 


(6.14) 


dx 
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which  has  second  moment  bounded  by 


MU', 


f  f*  -  M'"' 

Jo 


(«) 


dx  <  nw 


f<  „  u' j-) 

li,  •^■TwRTT 

ds  dx 


dx 


where  Doob’s  inequality,  Holder’s  inequality,  and  (6.8)  are  used.  Therefore 

From  (6.16)  with  u;„  =  1, 

EWhw^  =  o(i)[ /y‘Ei/,rd.sd.T 

L  JoJo 


Next, 


(6.15) 


(6.16) 


(6.17) 


(6.18) 


(6.19) 


=  B\!,fdsdx, 

by  Lemma  6.  Using  (6.8)  once  more, 

EWUW^  -  o{l}  JJ^  E\f^\-dsdx. 

It  can  be  checked  that 

I/ll  <  0{c^)\li  -  h\I{h  ^  0)  +  0(l)/(/i  =  0), 

V  *  w 

uniformly  in  t,  z.  Thus,  since  h  =  h  \i  h  ^  0, 

£^|/.|'  <  0{cl)E\h  -  /ip  +  0[\)P[h  =  0), 
so,  from  Lemmas  4  and  5, 

ff  E\h\Usdx  =  0{cln-^]. 

JoJo 

Combining  the  bounds  (6.16-6.20),  we  find  that  the  second  moment  of  the  Ihs  of  (6.10)  is  of 
order  (9(log<i„)^c^n“3^,  which  tends  to  zero  by  Condition  3.1.  This  establishes  (6.10). 

Step  2  We  now  prove  (6.11).  Let 

/t"'  (s,  u)/t~5(.<!,  t)  /?“*(s,  u)/i“2  (s.  .t) 


(6.20) 


6{s,u,x)  = 


j^h-Hs,v)dv  J^lr^s,v)dv 
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By  the  arguments  of  Step  1,  the  second  moment  of  ‘he  Ihs  of  (G.ll)  is  hounded  Ijy 

2  11  ^ 
j  Esupjj^  6(s,u,i)dA'(s,u)  dx  <  O(logd„)^^  \joL 

11  1  ■ 
<  o{hgd.)y’ii  E\S(s.  ii,x)\^ds  du  f/.rj 


=  O(logd„)^4« 

where  the  bound  on  the  triple  integral  is  from  Lemma  7.  □ 


0, 


Proof  of  Theorem  3.3  Define 


A’{t,z)  =  2f  f  X{s,x)f{s,x)dxds, 

Jo  Jo 

to  which  A{t,2)  converges  in  probability  under  the  general  alternative.  From  the  definition 
of  J  in  (3.7),  it  is  easily  checked  that 


dU{A-A-) 

didz 


h-^{t,z) 


A(t,z) 


y  lr'{t,u)X{t,u)du 
Iz  h-'^{t,v)dv 


Suppose  that  J(A  —  A*)  =  0.  Then  the  expression  inside  the  square  brackets  above  vanishes, 
so  that  ^  ^ 

X{t,z)  h~^{t,v)dv  =  h~\t,u)X{t,u)du, 

Taking  partial  derivatives  wrt  z  both  sides  gives  that  dX{t,  z)ldz  =  0  for  (<,  c)  €  [0, 1]  x  [O,/?], 
so  that  Ho  holds,  contrary  to  the  premise  of  the  theorem.  Thus,  J{A  —  .4'))  ^  0.  From 
arguments  in  the  proof  of  Theorem  3.2,  it  can  be  seen  that  ||(J  —  J){A  —  <4")|| — >0.  Hence, 


II  -  ^■)||-^||  J(>1  -  ^■)||  >  0.  (G.22) 

Along  the  lines  of  the  proof  of  Theorem  3.2,  it  can  be  shown  that  J{\/n{A  -  >4*))  converges 
weakly  in  Z)2([0)ll  (Oi/’D)  although  not  necessarily  to  Brownian  sheet,  cf.  the  proof  of 

Proposition  4.3  of  MU.  Similarly,  using  MU  (Proposition  3.2),  it  can  be  shown  ihat  J{y/n{A— 
A))  converges  in  the  same  sense.  The  triangle  inequality  gives 

-  ^))ll  +  -S'  +  PiMA  -  .4'))||  =  5  +  0/.(l), 

p 

the  equality  holding  by  the  continuous  mapping  theorem,  so  that  S — >co  by  (6.22).  Thus 
the  test  is  consistent.  □ 


Appendix 

The  following  lemma  is  routine. 

Lemma  1  Let  hjii,h2  be  Junctions  on  [0,  that  have  bounded  variation  and  arc  Lipschitz, 
with  h  nonnegative  and  bounded  away  from  zero.  Then  l/h,  \/Ti\  and  h\h-2  have  bounded 
variation  and  are  Lipschitz. 
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The  next  two  lemmas  collect  some  properties  of  weak  integrals  in  the  plane.  The  first  is 
a  version  of  the  integration  by  parts  formula.  Let  4> J  where  J'=  [o.a.']  x  [/».//]. 

Lemma  2  If  (  :  J  -*  R  is  cotitinuous  and  and  0(‘,b)  hurt  bounded  vurintiou. 

then  the  weak  integral  ffj<i>d(,  exists  and  is  equal  to: 

JJ^d{<l>{s,x)^(s,x))  +  JJ^^({s,x)d<i>(s,j  )  -  i{s,b')d<i>{s,b') 

+  /  ((s,b)d<f>(s,b)  -  f  ({a\x)d<f>{a.T)+  f  (,(a,x)  do{a,  x).  (.A.l) 

J  U  J  h 

Proof  Theorem  9.3  of  Hildebrandt  (1963)  gives  that  the  weak  integral 

JJ^  (4(^1  a;)  -  4(s,  b)  -  ^(a,x)  +  ^(a,  b))  d<^(s,  x)  (.4.2) 

exists,  and  coincides  with  the  weak  integral 

jj^{4>{s,x)  -  (f>{s,b')  -  4>{a\x)  +  d>[a\b'))di[s,x),  (A.3) 

which  exists  by  Theorem  8.8  of  Hildebrandt.  Theorem  5.8  of  Hildebrandt  shows  that  p(-,x) 
and  ^(s,  •)  are  of  bounded  variation  for  fixed  s,z.  (A.l)  can  then  be  obtained  by  rearranging 
the  terms  in  (A. 2)  and  (A.3).  □ 


Lemma  3  Let  (  be  a  stochastic  process  on  J .  If  the  weak  integral  ffj(i>df  exists  n.s.,  and 
the  stochastic  integral  J[j<i>d(  exists  in  the  L^-sense,  then  they  coincide  a.s. 

Proof  The  result  follows  immediately  from  the  definitions  of  the  stochastic  integral  and 
the  weak  integral,  and  the  fact  that  an  L^-limit  agrees  almost  surely  with  an  a.s. -limit.  □ 

The  next  lemma  is  a  refined  version  of  Proposition  3.3  of  MU,  giving  a  rate  of  convergence 
of  h  to  h. 


Lemma  4  There  exists  C  >  0  such  that 

/Y‘  E\h{t,  z)  -  h{t,  z)f  dt  dz  =  0{n-^ ). 

JoJo 

Proof  We  shall  use  much  of  the  notation  of  MU  (proof  of  Proposition  3.3),  without  re¬ 
defining  it  here.  As  in  MU, 

\h  -  hf  <  0{l)[\h  -  +\h-  -f  \h°  -  +  |/it  _  h‘f  +  (A..1) 

For  the  first  term. 


sup  \h{i,z)-h(t,z)f  < 
t,z 


(A. 5) 
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Also, 


Ar<">(<fs,x)  ] 

.  3 


r'/'  Ai"i(s,i)  ,  ,  i’ ,  . ,3,1  rr  m'-H'I.-.t) 


^  j:j:  "  H  -  j:j:  ^  - 


.  (A.G) 


The  first  term  in  (A. 6)  can  be  shown  to  be  bounded  using  (G.8).  The  expectation  in  the 
second  term  in  (A.G)  is  bounded  by 


(A.7) 


where 


^  /o  (V»(s,t))2 


and  we  have  suppressed  the  dependence  of  Mi{t)  on  x.  Let  [Mi]  and  (A/i)  be  the  quadratic 
variation  and  the  predictable  quadratic  variation  of  martingale  A/j,  respectively.  VVe  shall 
use  the  Burkholder-Davis-Gundy  inequality  (Dellacherie  and  Meyer,  1982,  p.2S7) 


E  sup  Mfiv)  <  CE[Mi]l 

u€(0,tl 


(A.8) 


Since  the  square  integrable  martingale  Mi  is  of  integrable  variation,  it  has  no  continuous 
part.  Hence  [A/i]<  =  (AAfi(u))^.  The  process 


at)  =  [Mi]t-{Mi),=  f 

Jo 


is  a  martingale,  so  E[Mi]l  <  2E{M\  )\  +  2£^^^(1)  =  2E{Mi)'\  +  2E{()^.  Ilut.  by  (6.8), 

It  then  follows  from  (A.7)  and  (A.8)  that  the  second  term  in  (.A.. 6)  is  of  order 

6- 

0{nWnfo{ - )  '  ->0, 

\nWnJ 

so  £",^^(1, 1)  <  oo,  and  from  (A. 5), 


£supi/t(^,c)-/t(<,x)r  = 

l.z 


(A. 9) 
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Since  sup  |/»(<,c)  -  h°{t,z)\  =  0(6„),  by  the  Lipscliitz  condition  on  /), 

(,I>(>n 


n'\h{t,2)-h°{t,z)\^dtdz=o{h„). 

JoJo 


(A.IO) 


For  the  third  term  in  (A. 4), 
sup  \h°{t,z)  -  h^{t,z)\ 


t.z 


= li^  /o‘  (^)  [/.'  (V)  s  '  ’ 


d. 


h{s,Xr)dx 


dx 
+  0(u;n) 


1 


<  -^0{wnb~^)  +  0(u;„)  =  0{Wnb~^). 

bn 

For  the  fourth  term  in  (A.4), 
sup 


(A.ll) 


t.2 

=  sup  E 
1.* 


^1  /■»  .,n-s 


Sl/o  "'(x 


'"'’■"'"“Xw-)* 


-Aik) 'Tf  At}lo'''(kr)k 


,Xr)  -  nWn 


{Yr{s)r 

Al")(s) 


3 

3 


/  1  \  ^  I  A 

<  of  -prr  )  {hndn)^  sup  E\h{s,Xr)  - 

\b‘^dr,J  S,r  I  (yv 


^  ds 


AW(s) 


Using  the  Lipscliitz  property  of  A  and  (6.8),  the  supremuni  term  above  is  bounded  by 


0(u;^)  +  0(1)  sup  £ 


n(")(s) 


<0(u;^)  +  0(l) 


0 


nw„ 

1 


-  f{s,Xr) 


\nWn 


+  0(fc-^"'"" )  =  0 


+  0(i)supP(y;f")(s)  =  0) 

s.r 

1 


liVK. 
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where  we  Imve  used  the  Lipschitz  property  of  /  and  llie  fact  tliat  }'/"*(.•<)  is  a  i)inoniial  r.v. 
with  mean  of  order  0(niu„)  to  bound  its  fourth  central  moment.  Tliereforc 

=  0  (1)  0  0  ( (A.li) 

Finally,  for  the  fifth  term  in  (A. 4), 

sup  E\R(Uz)\^  <  sup[ER\t,z)]l  (A. 13) 

(,Z  1,2 


nn  .  -  ”  ^  r  >.Y<  -  dA/<’‘)(.s) 

6„  j/o  (  b„  / 

Then  R{t,z,-)  is  a  n  rtingale.  Using  the  same  arguments  that  were  applied  to  .\/i, 

ER*{t,z)<E  sup  /?"(t,s,u)  <CE[/?(t, ->,•)]?<  0(l)(£’(/?(<, --,•))?  + (A. 14) 

0<u<I 

where  ^(u)  =  [R{t,  z,  •)]„  —  r,  •))^  is  a  martingale.  Since  no  two  of  the  counting  processes 
r  =  1, 2,  •  •  • ,  d„  jump  simultaneously, 

v<u 

^  /  n  ■2fz  ~Xr\, ,2  ft  —  v\AN^’'^{v) 


-  yyfJL.)  (LzJLi izl ^ 


f  ^  2*2/*  /**  i’2f^ 

\wJ  S"'  mi. 


It  follows  that 


The  last  two  terms  in  (-4.14)  can  be  bounded  above  as  follows: 


-  °[wJ .S mi 

=  o(i)(.A;)-'. 


(A. 1C) 
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supE^^ll)  =  sup£;(0, 

t,z  C,z 


=  0(Vi(nHl)-'. 

From  (A.13-A.17),  we  get 

sup£|fl((,.-)|“<0(l)k«r‘. 

<.2 

Finally  combining  the  bounds  for  the  five  terms  in  (A. 4),  we  have 

//  E\h{t,2)-  h{t,2)f  dtdz  =  0(u)^6;®  +  6n  +  ^i’n*;®  +  (ni^nin)”^  +  ('^^n)' 

Jo  Jo 

=  0{n-<) 

by  Condition  3.1,  where  C  =  min(/?,y(a  —  ^))  >0.  □ 


(A.17) 


(A.IS) 


Lemma  5 


P{h{t,2)=^0)dld2  =  0{n-<). 

JoJo 


Proof  Let  0<c<C<oobe  lower  and  upper  bounds  for  h. 
P{h{t,z)  =  0)  =  P{h{i,2)  <  c~^)  + P{h{t,2)  >  Cn) 


<  P{h{t,  z)  -  h{t,2)  <  C-'  -  c)  +  P{h{t,  2)  -  h{t,  2)  >  c„  -  C) 
=  P{h{t,z)  -  h{t,z)  >c-  c"’)  +  PCh{t,2)  -  li{t,2)  >  c„  -  C) 

^  E\h{Uz)-h{t,z)f  E\h{t,2)-h{t,2)f 

(c-c;‘)3  (c„-C)3 


so  the  result  follows  by  Lemma  4.  □ 


Lemma  6 


Proof  The  expression  on  the  Ihs  above  is  bounded  by 


supi;[s^OKK(>'‘"’(^-)  7^  0)  +  c>(v^);(r‘">(Lz)  =  o) 

<  0{Jnwl)  +  0(77)sup  P(Y^''\i,  ')  =  0) 

t,: 


where  the  last  inequality  comes  from  (6.9).  □ 
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Lemma  7  For  S{s,u,x)  defined  by  (6.21), 


s,u,u:)}i^dsdudx  =  0{cl,^ 


Proof  Let 


7i  =  l(k{s,u)  ^  0,  li{s,x)  j  ¥"  0): 

h  =  l(k{s,  u)  =  G)  +  l{h{s,x)  =  0)  +  7(  J'  r’(.s  v)  dv  =  0). 


Then 


|5(s,u,x)|  =0(Cn) 


(h~\s,u)  -  h~\s,u))h~^{s,x)  J  h~\s,v)dv 
+(h~^{s,x)  —  h~^(s,x))h~^{s,u)  J  h~^{s,v)dv 
+h~^{s,x)h~^{s,u)  J  {h~^{s,v)  —  h~^{s,v))dv 


1\  +  0(12) 


=  0(c„) 


c,^,lh~'(s,u)  -  h  ‘(s,u)l  +  |/rv(.s.r)  -  /r^.s,ar)| 


+  J  \h~^{s,v) -h~^{s,v)\dv 


A +0(72) 


=  0(Cn)  |c| \k{s,  u)  -  h{s,  u)l  +  Cn\h{s,  x)  -  h{s,  x)| 

+c„  J  \'h{s,v)  -  h{s,v)\I{h{s,v)  0)dv 


+ 


J  l{h{s,  v)  =  0)  dv 


+  0(A). 


V  A  w 

where  we  have  used  the  fact  that  h  =  h  when  h  ^0.  Thus, 

JS|6(s,U,x)P  <  0(c„)^^Cn£^|/l(s,u)  -  /j(s,u)p  +  C^£|/t(s,x)  -  /l(s,x)p 

+c^  /  E\h{s,v)  —  h{3,v)\^ dv  +  f  P{h{3,v)  =  0)  dv 
Jo  Jo 

+0(1)[p(M5,u)  =  0)  +  P{h{s,x)  =  0)  +  p(^J^  lr^{s,v)dv  =  0)  . 

Notice  that  if  h{s,  1)^0,  then  c”*  <  li{s,  1)  <  c„.  Since  h{s,  •)  is  continuous,  c"*  <  h{s,  v)  < 
c„  in  a  small  interval  v  £  [1  —  1],  £  >  0.  Hence  c"-  <  h~^{s,  v)  <  c„,  for  v  €[!—£,  1], 

Hence  /j  h~^{s,v)dv  ^  0,  for  x  <  1.  Therefore, 


■(/: 


h~^(s,v)  dv  =  0]  <  P{h{s,  1 )  =  0),  for  x  <  1. 


Now  apply  Lemmas  4,  5  and  the  fact  that  /J  P{h  =  0)dt  =  0{n  *•),  which  is  easy  to 

see  from  proofs  of  Lemmas  4  and  5,  to  complete  the  proof. 
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